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Abstract: We consider extensions of the Standard Model with a number of additional 
heavy quarks which couple to the Higgs boson via top-like Yukawa interactions. We con- 
struct an effective theory valid for a Higgs boson mass which is lighter than twice the lightest 
heavy quark mass and compute the corresponding Wilson coefficient through NNLO. We 
present numerical results for the gluon fusion cross-section at the Tevatron for an extension 
of the Standard Model with a fourth generation of heavy quarks. The gluon fusion cross- 
section is enhanced by a factor of roughly 9 with respect to the Standard Model value. 
Tevatron experimental data can place stringent exclusion limits for the Higgs mass in this 
model. 



1. Introduction 



The discovery of the Higgs boson will introduce a new era in particle physics. The long 
standing theoretical problem of understanding the mechanism for electroweak symmetry 
breaking will be tackled for the first time using direct experimental findings. The measure- 
ment of the Higgs boson mass and the production cross-sections of its various signatures will 
be important constraints in formulating a theory of particle interactions at high energies. 

The interaction of the Higgs boson and gluons is particularly important at hadron 
collider experiments. In the Standard Model (SM), the gluon fusion cross-section is the 
largest among all production cross-sections. The LHC will be able to discover the SM 
Higgs boson in this production channel for the full range of its allowed mass values. The 
branching ratios for the decays of the SM Higgs boson are dominated by other Higgs boson 
interactions which involve the bottom quark and electroweak gauge bosons. However, the 
Higgs-gluon interaction is still not negligible; a fraction of up to about 7% of Higgs bosons 
may decay to gluons, depending on the Higgs boson mass. 

Given that gluons are massless, a Higgs-gluon interaction arises as a loop effect via 
other massive coloured particles which couple to the Higgs boson. Physics beyond the 
Standard Model can alter significantly the strength of this interaction in various ways. 
One possibility is that new coloured particles are not much heavier than the top quark, 
and their contribution is therefore not suppressed. A second possibility is that new coloured 
particles may be heavier than the top-quark but they have an enhanced Yukawa coupling 
to the Higgs. A third possibility is that such particles are quite heavier than the top-quark, 
but their multiplicity is large, thus building up a significant cumulative contribution. 

A Higgs boson is often assumed to be lighter than about twice the mass of the top-quark 
and twice the mass of new undiscovered particles which are hypothesized in extensions of 
the Standard Model. Light new particles are hard to accommodate given the vigorous 
experimental searches for new physics at LEP and the Tevatron. On the other hand, they 
cannot be very heavy or, alternatively, they must have a rather strong Higgs coupling if 
they contribute in reducing the fine tuning of the Higgs mass. Therefore, it is important 
to calculate their contribution in the gluon fusion process as well as in the decay of a Higgs 
boson to gluons. 

Assuming a Higgs boson which is lighter than production thresholds of new heavy 
particles, we can factorize the effect of new physics and QCD in the process gg H, 
by means of an effective field theory where the top-quark and all other possible heavy 
coloured states which couple to the Higgs boson are integrated out. The effect of these 
heavy particles is included in the Wilson coefficients of an effective theory with operators 
of the Higgs boson and light quarks and gluons. 

The possibilities for viable extensions of the Standard Model which alter the Higgs- 
gluon interaction are many, and an equal number of matching calculations is required for 
their study. This is a rather easy task at leading order in the strong coupling. However, 
experience from the Standard Model shows that a precise estimate of the gluon-fusion 
cross-section and the Higgs decay width to gluons requires calculations through next-to- 
next-to-leading-order (NNLO) in the strong coupling. 
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In this paper, we consider extensions of the Standard Model with additional heavy 
quarks. We assume that these quarks have a Higgs Yukawa interaction of the same type 
as Standard Model quarks. The existence of such quarks has dramatic implications for the 
Higgs production cross-section in gluon fusion. At leading order, and in the limit where 
the heavy quarks are much heavier than half the mass of the Higgs boson, the cross-section 
scales as n|, where is the number of heavy quarks. Current measurements at the 
Tevatron [1] and early data from the LHC can therefore constrain severely such models. 
We first construct an effective Lagrangian integrating out the top-quark and the additional 
heavy quarks. We compute the Wilson coefficient of the Higgs-gluon effective interaction 
through NNLO in the strong coupling expansion. Finally, we present numerical results 
for the gluon fusion cross-section at the Tevatron in a specific model with a fourth quark 
generation. 



2. The effective Lagrangian 

We consider an arbitrary extension of the SM throTigh new heavy qiiarks transforming 
under the fundamental representation of the QCD gauge group SU{3). The number of 
heavy quarks, including the top, is nh- We will denote their mass by m^, with q = 1 . . . n^- 
We assume that the new quarks, as the SM top, couple to the Higgs boson H through their 
mass. Therefore, the Lagrangian we begin with is 

^ = ^QCD + ^%{'i'P-^q)%+ , Cy = ^mqtjjglpq. (2.1) 

q=l q=l 

Here is the covariant derivative in the fundamental representation and CJqcd 
QCD Lagrangian with only the ni flavours of light quarks. We take these quarks to be 
massless. 

We focus on the changes that the heavy quarks induce on the Higgs production through 
gluon fusion. When the quarks that couple to the Higgs boson are heavier than half the 
Higgs boson mass, we can integrate them out. In this limit, we can replace the original 
Lagrangian (2.1) with an effective Lagrangian 

= C^Sd -Ci^O,. (2.2) 

Ci is the Wilson coefficient [2] relative to the only dimension-four local operator Oi that 
arises when we integrate out the heavy quarks and all the quarks remaining are massless [3] , 

= i G'l^G""''' . (2.3) 

In this expression, G'^j, is the field strength tensor in the effective theory. In Eq. (2.2), 
^QCD describes the interactions among light quarks. It has the same form as jOq^^, but 
with different parameters and field normalizations because of the contributions from heavy 
quarks loops. We relate the parameters in the effective theory to the parameters in the 
full theory through multiplicative decoupling constants Q. We will denote quantities in 



- 2 - 



the effective theory with a prime. The derivation of the decouphng constants is reviewed 
in [4]. In Section 4, we describe the main steps of their calculation and give the relevant 
results. 



3. Method 

We compute the Wilson coefficient Ci up to three loops. Diagrams containing both the 
heavy mass scales appear for the first time at the three-loop order. We start from the bare 
amplitude M^g^jj for the process gg ^ H in the full theory, 

M'gg^H ^ M'^X'(P^^P^)CA^ ■ (3-1) 

Here, pi and p2 are the momenta of the two gluons with polarizations Cal and £02- This 
amplitude is related to the bare Wilson coefficient by [4] 

- (7V2-l)(d-2)(pi.p2)2 ^^^^^ iPi,P2\,=,,=o ■ (3-2) 



N is the number of colours and d = 4 — 2e is the dimension of space-time. Bare quantities 

53 



are denoted by the superscript "0". The factor ("3 is the bare decoupling coefficient by 



which the bare gluon field G'J^'" is rescaled in the effective theory. 



G'T = \/Q G'^^ . (3.3) 

We generate the Feynman diagrams T for the amplitude through three loops using 
QGRAF [5]. We then perform an expansion of all diagrams in the external momenta 
Pi,P2, by applying the following differential operator [6] to their integrand: 

00 

^•^ = J2^P1 ■ P2r Pn-^]p,=p2=0 ' (3-4) 



n=0 



with 

Vo = l, v,= '-n,2, ^2 = -^^^3T)^{°nn2.-ciDM, (3-5) 
and Dij = g^'^Q^- 

Differential operators of higher orders are not needed for the expansion in the external 
momenta at leading order. 
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After Taylor expansion, all the Feynman diagrams are expressed in terms of one-, two- 
and three- loop vacuum bubbles by using linear transformations on the loop- momenta ki'. 



with 



^iN=y-^^' (3-6) 

JSa [i^l, 1^2,1^3, 1^5, 1^6, l^7\ = J ^.^d/2)3 V^^V^V^V^V^^ ' ''^"^^ 

, r i_ f d^kid%d% 1 

J3& [i^i, 1^2, 1^3, 1^4, i^5, z^ej = J — ^—d/rjs — v^^v^vfvfV^V^ ' ^ ' 

/3c[.l,.2,.3,.4,.5,.6] = J ^^^,^,^3 ^ . ^ ^ .3 ^ p p ' (^"lO) 



= A^f - m2 , 

= fci - m2 , ^2 = fei - m2, , 

= fci - , P3 = kl- ml, , 

= {ki -k2 + k^f - ml , (3.11) 

= (fci - k2)^ , 

n = {h - k^f , 

V7 = {k3 - k^f , 



and fj, Vi positive or negative integers. The third three-loop vacuum bubble I^c contains 
two heavy quarks of different mass, mq and m^/. 

We perform a reduction of the above integral topologies to master integrals using the 
algorithm of Laporta [7] and the program AIR [8]. We find five master integrals, 

1i = h[l\ 

= -K)^-^r(-i + 6), (3.12) 

X2 = /3a[l, 0,1, 1,1,0] 

_ / 2x2-3e r^(l - 6)r(6)r^(-l + 2e)T{-2 + 3e) 
" ^'^'i> r(2-e)r(-2 + 4e) ' 

X3 = /36[l, 1,1, 1,0,0] , (3.14) 

X4 = /3c[l, 1,1, 1,0,0] , (3.15) 

X5 = /3c[2, 1,1, 1,0,0] . (3.16) 



Single-scale master integrals appear in the calculation of the SM Wilson coefficient [9, 10] 
and can be computed with MATAD [11]. For the remaining two-scale master integrals we 
used result from [12]. We checked all the master integrals independently through sector 
decomposition with the program of Ref. [13]. 
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After these steps, the RHS of Eq. (3.2) becomes 



9=1 



1 ra% 



TT 



-2e' 



-1 + e [l + 21ogl 



log2(m°) + log(m°) + 



TT 



+ 



+ 



IT 
TT 



log(m°) + 



31 
36 



24 



1 1 



31og(m°) 223 



16 



576 



+ rii 



5 1 103 51og(m°) 
144 e 864 24 



-— log2(m°) + 



2231og(m°) TT^ 5975 



16 



96 



\J2 



TT 



1 1 

16? ~ 16e 



128 ^ 3456 

89' 



3(log(m°) + log(m°,)) + 



+ 0{e 
1 



18 



log(m°)log(m°,) 



+^ NK) + i°g«)] + ^ + i°g'«)] + ^^^^^^ + 



16 



1728 



(3.17) 



The first sum in this expression runs over all single-scale diagrams and corresponds to Uh 
copies of the SM Wilson coefficient. The second sum accounts for the 3-loop diagrams in 
which either of the two massive quarks couples to the Higgs boson. We already symmetrized 
it over q,q'. In Eq. (3.17) we introduced the factor 



Se = e 



(3.18) 



4. Decoupling and renormalization 

The RHS of Eq. (3.17) contains the bare masses of the heavy quarks and the bare 
coupling constant in the full theory; Ci = Ci{a^,mg). The bare strong coupling in 
the full theory is related to the bare strong coupling in the effective theory by the 
decoupling constants (g [4, 10] , 

«;°=(c°)^«° . (4.1) 

Similiarly, 

a, = {Cg?as . (4.2) 

Using these relations, we obtain the bare Wilson coefficient as a function of the bare 
parameters in the effective theory and of the bare mass of the heavy quarks in the full 
theory, Cj* = C'l{a'^ ,171^). The bare parameters are related to the renormalized ones 



-5- 



through multiplicative renormalization constants Zi as 



(4.3) 
(4.4) 



All the parameters in Eq. (4.3) arc in the effective theory and all the parameters in Eq. (4.4) 
are in the full theory. Finally, we renormalize the Wilson coefficient itself through a renor- 
malization factor Zii [3,14,15], 



Ci 



1 



(4.5) 



4.1 Details of the calculation 

A convenient way to compute the gluon field decoupling (^3 is through the relation 

Cl = l + U%{p = 0), (4.6) 

where Hq is the transverse component of the gluon self-energy in the full theory. This 
quantity is computed at zero external momentum. Since we work in dimensional regu- 
larization, only diagrams containing at least one massive quark contribute. We find only 
one diagram per heavy flavour at one loop and seven at two loops. We employ the same 
calculation techniques as in Section 3. Our result reads 

1 log(m°) 7r2 + 241og2(m°) 




6e 



+ e- 



72 



1 



91 



l + 241og(m°) 3 o n X 

+ - Iog2(m0) + — log(mO) + 

64e 4 ^ ^ 16 1152 



+ 



TT 

64 



(4.7) 

The bare decoupling parameter of the strong coupling constant, can be computed 



as 



So 



(4.8) 



where (f, Q and are the bare decoupling constants of the gluon-ghost vertex, of the 
ghost field and of the gluon field respectively. 

The bare decoupling constant of the ghost field is calculated from the ghost self-energy in 
a similar way as ^3. At one loop, there is no diagram contributing to the ghost decoupling. 
At two loops, there is only one diagram per heavy flavour. We find 



TT 



2 n,j 

E 

q=l 



3 1/5 3 , , o^ 

n + \ log(m^) 

64e2 e V 128 16 ^ 



768 



-^logK)-^log2(m°) 



(4.9) 
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The decoupling of the gluon-ghost vertex (^^ is given by 

CT = i + r?Gr,(o,o), (4.10) 

where T^Q^{p,p') is extracted from the IPI amputated gluon-ghost Green function and p 
and p' are the incoming four-momenta of fj and G respectively. Again, this term receives 
no contribution at one loop. At two loops, there are 5 non-masslcss diagrams for each 
heavy flavour. Two of them vanish because of colour, and the other three add up to zero. 
Therefore 

CT = 1 + 0{al) . (4.11) 
Inserting Eqs. (4.11, 4.9, 4.7) into the relation (4.8) we find 



c: 



1+ 



+ 



IT 



TT 



vie 6 



96e2 



24e 



6 



144 



24 



24 



Here we introduced the notation 



9=1 



L°, = Elog2(m°) 

9=1 



(4.12) 



(4.13) 



We now renormalize the mass of the heavy quarks in the full theory according to 
Eq. (4.4). The mass renormalization constants in the full theory are obtained from 
the one- and two-loop corrections to the quark propagator [16]. We review here the main 
steps of this calcTilation. Let us denote the sum of all the one-particle irreducible (IPI) 
insertions into the quark propagator as — zSo(p), 



-zSo(p) = -zSj^(p)-zE^^(p) + ... , 



ILi 



(4.14) 



where —iY^ij)) in the sum of all the one-loop IPI diagrams in the bare theory and so on. 
The full quark propagator then reads 



mi 



We can split So(p) as 



Mp) = 5^10 (/) + (^- m°)S2o(/) ; 



conversely, the quantities ^Sio, S20 are extracted from the bare self-energy Eq as 



mo 



— Eio = -'IV — So + ^So 



mo 



^20 



V^(^^^o) 



Combining Eqs. (4.15) and (4.16) we obtain 



m° = mqi/i) 



1 - — S}^ - I —Ef'^ + — Si'^S 



mn 



2L 



mn 



mn 



(4.15) 
(4.16) 

(4.17) 
(4.18) 



-7- 



The RHS of Eq. (4.18) is expressed in terms of renormalized quantities. On the other 
hand, the one-loop result for Y1\q originally depends on the bare gauge parameter and on 
the bare strong coupling and is divided itself by the bare mass. Therefore, one also needs 
the one-loop renormalization of these three quantities in order to compute through 
two loops. 

We compute the self-energy for arbitrary external momcntTim p, and set the mass of the 
heavy quarks to zero. In the MS scheme we recover the result [16] ^ 

2r. 2n/ 



^m„ — J- r 

TT e 



TT 



1 /45 



24 



1 



101 



+ -i-l^+72"^ 



(4.19) 



This relation holds both in the full theory, where the number of active flavours is 
Hf = ni + nh, and in the effective theory. In this case, Uf = ni and we need to replace the 
renormalized strong coupling in the full theory with the one of the effective theory. 
So far, the coefficient for the mass decoupling and renormalization still depends on the 
renormalized strong coupling in the full theory. We decouple it using the relation [4] 



-a.s{\i) = dasifi) ■ 



(4.20) 



The strong coupling renormalization constants are related to the coefficients of the /3 func- 
tion as 

2 



+ 



TT 



TT 




(4.21) 



Here /?q and denote the first two coefficients of the /? function in the light-flavours theory. 



11 



16 



38 

102 m 

3 



(4.22) 



Combining Eqs. (4.4 -4.21) we get 



m = mq{n) i 1 



with 



TT e 



TT 



Lr. 



Ab-2{nh + ni) , 48L^ + 10(n/, + n/) - 303 



q=l 



24e2 



+ 



144e 




The next step is the renormalization of the bare strong coupling in the effective theory. 



ag. The relevant renormalization constant is given in Eq. (4.21). 

We finally renormalize the bare Wilson coefficient Ci{a'g,mq) using [3,14,15] 

1 



1 + 



da's (a*) 



log ^4 



TT 



/3( 



(4.25) 



^Note that in our conventions d = 4 — 2£, while in Ref. [16] d = A + 2e. This explains the sign difference 
in the 1/e terms. 
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Our final result for the renormalized Wilson coefficient reads 

2 



c 



3 



TT 



TT 



TT 



1877 
l92 



■rih 



67 



77 2 , 19L^ 
576"'^ + ^ 



ILr. 



(4.26) 



This is the main result of our paper. Note that the second term in the square brackets 
comes from diagrams containing two massive quarks loops. For n/j = 1 we recover the SM 
Wilson coefficient [9,10,17,18] through order 0{a'^). 

At leading order in the heavy-quark expansion, and assuming a massless bottom quark, 
the gluon fusion cross-section and the decay width of the Higgs boson to gluons are pro- 
portional to the square of the Wilson coefficient. In this limit, their ratio with the corre- 
sponding quantities in the Standard Model are: 



gg)inH) 



a{gg H)iSM) y{H gg)^^^) 



TT 



rih 



—n in 



1) + 



4 19\v-n 

-n, + -J)^log 



rutin) 



(4.27) 



where mt the mass of the top-quark. The 0{ag) term in the above expression is generally 
very small. However, in a realistic phenomenological study [19, 20] the exact quark mass 
dependece of the cross-section as well as effects due to electroweak corrections need to be 
accounted for through NLO. 



5. Numerical Results for gluon fusion cross section at the Tevatron 

In this Section, we present our numerical results for the cross-section at the Tevatron, in 
a Standard Model with four generations. The Wilson coefficient for the fourth generation 
model is obtained by considering three heavy quarks in Eq. (4.26). We set the top-quark 
mass to 

mt = 170.9 GeV . 

For the fourth generation we consider two scenarios, corresponding to fourth generation 
down-quark masses of: 

ruB = 300 GeV, = 400 GeV 



and an up-quark mass given by 



mT-mB = 50 GeV + 10 log GeV. 

V 115 GeV/ 



(5.1) 



This choice is permitted by constraints from electroweak precision tests, as described in 
Ref. [21]. 

The calculation of the total cross-section differs from the Standard Model only in the 
expression of the Wilson coefficient for the low energy effective Lagrangian. We combine 
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Eq. (4.26) with the known results for the Standard Model total cross-section at NNLO in 
the large top- mass limit of Rcfs [22-24] . 

Adopting the same approach as in Ref. [19], we first compute the ratio of the NNLO 
and LO cross-section in the effective theory. We estimate the contribution from Feynman 
diagrams with only top-quark and fourth generation quarks to the total cross-section, by 
multiplying this ratio with the exact leading order contributions of heavy quark diagrams 
in the full theory. 

NNLO;{t,B,T) _ / ^ \ LO;{t,B,T) ,^ 

"heavy — I ^LO;{t,B,T) \ "exact W-^^ 

\ / effective 

These contributions are enhanced by roughly a factor of 9 with respect to the corresponding 
Standard Model results, since 

LO;{t,B,T) _ LO;{t) . ^ 

"exact — ^ "exact W-"^/ 

within a few percent. 

Contributions from diagrams with bottom quark loops are small and we compute them 
exactly through the NLO order approximation [25,26]. In comparison to their Standard 
Model counterparts, the most important interference terms of diagrams with bottom quarks 
only and diagrams with any of the heavier quarks are only enhanced by roughly a factor of 
three. Therefore, these contributions are suppressed by roughly a factor of ~ 3/9 in this 
model. 

Finally, we include two-loop electroweak corrections from light-quark loops from the 
first two generations [27] in the complex mass scheme, and the corresponding three-loop 
mixed QCD and electroweak corrections as in Ref. [19]. These are enhanced by a factor of 
roughly 3 in comparisoin to the Standard Model, and are therefore suppressed by a factor 
of 1/3 in this model. Wc ignore electroweak corrections with quarks from the third and 
fourth generation, which arc already found to be very small in the Standard Model [28] for 
the Higgs mass range accessible at the Tevatron, when a complex mass scheme is employed. 

In Table 1, we present the cross-section at a renormalization and factorization scale 
IJ, = Hf = fir = mn/^ and estimate the scale variation error by varying the common 
scale /X in the interval [^,mi/]. We use the MSTW2008 NNLO parton distribution 
functions [29] , and compute the uncertainty (with 90% CL) to the cross-section due to 
the parton distribution functions (including the parametric uncertainty of the value of a^) 
according to Ref. [30]. 

The scale variation uncertainty and the uncertainty from the parton distributions are 
the dominant uncertainties. Essentially, they are the same as in the Standard Model cross- 
section with only three generations. We note that Ref. [19] preceded the release of Ref. [30] 
where it became possible to include the uncertainty of the as value in the fitted parton 
densities and the corresponding parton density uncertainty was estimated to be smaller. 

6. Conclusions 

In this paper, wc constructed an effective field theory for extensions of the Standard Model 
with many heavy quarks coupling to the Higgs boson via top-like Yukawa interactions. We 
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have computed the required Wilson coefficient of the —^Gfj^vG^'^ operator through NNLO 
in the strong couphng expansion. We found the result 



Ci = 



3 TT 4 TT 



TT 



1877 



192 576 

/67 2L„ 



77 2 , 



(6.1) 



This result can be utilized by the ongoing experimental studies at the Tevatron and the 
LHC to constrain such models. The cross-section and the decay width are enhanced by 
roughly the square of the number of heavy quarks with respect to the Standard Model. 
The Tevatron has put stringent limits on the Standard Model Higgs boson gluon fusion 
cross-section [1]. Equivalent studies can be performed in models with additional quarks. 
We have presented numerical results for the gluon fusion cross-section at the Tevatron in 
a non-minimal Standard Model with four generations. The theoretical uncertainties of the 
cross-section are practically independent of the number of heavy quarks and very similar 
to the Standard Model. 
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Table 1: The NNLO cross-section for Higgs production via gluon fusion at the TEVATRON. cr(i) 
corresponds to tub = SOOGeV and C7(2) to ms = 400GeV. The mass of the fourth generation up 
quark is given by Eq. (5.1) 
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